Rotating Stars in Relativity 



Nikolaos Stergioulas 



OO 



Department of Physics, University of Wisconsin- Milwaukee 
PO Box 413, Milwaukee, Wl 53201, USA 
and 

Max Planck Institute for Gravitational Physics 
(Albert-Einstein-Insitute) 
D-14473 Potsdam, Germany 
email: niksterg@aei-potsdam.mpg.de 



(N 
> 
(N 

O 

o 

00 
ON 

a"- 

I 

?H 



Review article for Living Reviews in Relativity 
Ettp : / /www . livingreviews . or^ 



Abstract. Rotating relativistic stars are receiving sig- 
nificant attention in recent years, because of the informa- 
tion they can yield about the equation of state of mat- 
ter at extremely high densities and because they are one 
of the more possible sources of detectable gravitational 
waves. We review the latest theoretical and numerical 
methods for modeling rotating relativistic stars, includ- 
ing stars with a strong magnetic field and hot proto- 
neutron stars. We also review nonaxisymmetric oscil- 
lations and instabilities in rotating stars and summarize 
the latest developments regarding the gravitational wave- 
driven (CFS) instability in both polar and axial quasi- 
normal modes. 



1 Introduction 

Rotating relativistic stars are of fundamental interest in 
physics. Their bulk properties restrict the proposed pos- 
sible equations of state for densities larger than nuclear 
density. Their oscillations can become unstable, produc- 
ing gravitational waves that could be detectable, provid- 
ing thereby a new way of probing the interior of neutron 
stars. 

Recent research has considerably advanced our under- 
standing of these objects. There now exist several inde- 
pendent numerical codes for obtaining accurate models of 
rotating neutron stars in full general relativity. Three of 
these codes have been shown to agree with each other to 
remarkable accuracy and one code is available as public 
domain for use by other researchers. 

The numerically constructed maximum mass models, 
for different proposed equations of state, differ by as much 
as a factor of two in mass, radius and angular velocity, a 
factor of five in central density and a factor of eight in the 
moment of inertia. These large uncertainties show that 



our understanding of the properties of matter at very high 
densities is currently rather poor. 

Despite the different maximum rotation rates, corre- 
sponding to different candidates for the equation of state 
of neutron-star matter, one can place an absolute up- 
per limit on the rotation of relativistic stars, by imposing 
causality as the only requirement on the equation of state. 
It then follows that gravitationally bound stars cannot ro- 
tate faster than 0.28 ms. 

Although observed magnetic fields in neutron stars 
have a negligible effect on neutron-star structure, a suf- 
ficiently strong magnetic field acts as a centrifugal force 
on a relativistic star, flattening its shape and increasing 
the maximum mass and rotation rate for a given equa- 
tion of state. The magnetic field strength of a stationary 
configuration has been shown to have an upper limit of 
B - 10" G. 

Rapidly rotating proto-ncutron stars arc shown to have 
an extended envelope, due to their high temperature and 
the presence of trapped neutrinos. If the equation of state 
is softened, as the neutron star cools, by a large ampli- 
tude phase transition, then the nascent neutron star may 
collapse to a black hole. A surprising result is that a 
supramassive proto-neutron star, even though it contracts 
during cooling, evolves to a cold neutron star of smaller 
angular velocity. 

In rotating stars, nonaxisymmetric perturbations have 
been studied in the Newtonian and post-Newtonian ap- 
proximations, in the slow-rotation limit and in the Cowl- 
ing approximation but fully relativistic quasi-normal 
modes (except for neutral modes) have yet to be obtained. 
The effect of rotation on the quasi-normal modes of os- 
cillation is to couple polar and axial modes and to shift 
their frequencies and damping times causing some modes 
to become unstable. 



Nonaxisymmetric instabilities in rotating stars can be 
driven by the emission of gravitational waves (CFS- 
instability) or by viscosity. The onset of the CFS- 
instability has now been computed for fully relativistic, 
rapidly rotating stars. Relativity has a strong influence 
on the onset of the instability, allowing it to occur for less 
rapidly rotating stars than was suggested by Newtonian 
computations. 

Contrary to what was previously thought, nascent neu- 
tron stars can be subject to the I — 2 bar mode CFS- 
instability, emitting strong gravitational waves. The fre- 
quency of the waves sweeps downward through the op- 
timal LIGO sensitivity window and first estimates show 
that it could be detectable out to the distance of 140 Mpc 
by the advanced LIGO detector. 

The viscosity-driven instability is not favored by gen- 
eral relativity but, as a new relativistic computation 
shows, is absent in rotating neutron stars, unless the equa- 
tion of state is unexpectedly stiff. 

Axial fluid modes in rotating stars (r-modes) received 
renewed attention since it was discovered that they 
are generically unstable to the emission of gravitational 
waves. The r-mode instability can slow down a newly- 
born rapidly rotating neutron star to Crab-like rotation 
rates. First results show that, the gravitational waves 
from the spin-down (directly, or as a stochastic back- 
ground) could be detectable by the advanced LIGO or 
VIRGO detectors. 

The present article aims at presenting a summary of 
theoretical and numerical methods that are used to de- 
scribe the equilibrium properties of rotating relativistic 
stars and their oscillations. It focuses on the most recently 
available preprints, in order to rapidly communicate new 
methods and results. At the end of some sections, the 
reader is pointed to papers that could not be presented 
in detail here. As new developments in the field occur, 
updated versions of this article will appear. 



2 The Equilibrium Structure of 
Rotating Relativistic Stars 

2.1 Assumptions 

Although a relativistic star has a complicated structure 
(solid crust, magnetic field, possible superfluid interior 
etc.), its bulk properties can be computed with reasonable 
accuracy by making several simplifying assumptions. 

The matter is modeled as a perfect fluid because obser- 
vations of pulsar glitches have shown that the departures 
from perfect fluid equilibrium due to the solid crust are of 
order 10~^ Q. The temperature of a cold neutron star is 
assumed to be K because its thermal energy (<< IMeV 
^ 10^" K) is much smaller than the Fermi energy of the 
interior (> 60 MeV). One can then use a zero-temperature 
(one-parameter) equation of state (EOS) to describe the 



matter: 



e = 6(P), 



(1) 



where e is the energy density and P is the pressure. At 
birth, a neutron star is differentially rotating but as the 
neutron star cools, shear viscosity, resulting from neutrino 
diffusion, aided by convective and turbulent motions and 
possibly by the winding-up of magnetic fleld lines, en- 
forces uniform rotation. At present, it is difficult to accu- 
rately compute the timescale in which uniform rotation is 
enforced, but it is estimated to be of the order of seconds 

Within roughly a year after its formation, the neutron 
star temperature becomes less than lO^K and its outer 
core becomes superfluid (see |^ and references therein). 
Rotation causes the superfluid neutrons to form an array 
of quantized vortices, with an intervortex spacing of 



d„ - 3.4 X lO^^f^, ^/^cm, 



(2) 



where is the angular velocity of the star in 10^s~^. On 
scales much larger than the intervortex spacing, e.g. on 
the order of 1 cm, the fluid motions can be averaged and 
the rotation can be considered uniform [Q.The error in 
computing the metric is of order 



10" 



-11 



(3) 



where R is a typical neutron star radius 

The above arguments show that the bulk properties of 
a rotating relativistic star can be modeled accurately by 
a uniformly rotating, zero-temperature perfect fluid. 

2.2 Geometry of Space-Time 

In relativity, the space-time geometry of a rotating star in 
equilibrium is described by a stationary and axisymmetric 
metric of the form 



ds 



e'^dt' + e'^id^ - udt f + e^"{dr' + r^d9^), (4) 



where v, ijj, uj and a are four metric functions which de- 
pend on the coordinates r and 9 only (unless otherwise 
noted, we assume c = G = 1). The perfect fluid has a 
stress-energy tensor 



T"'' = (e + P)u%'' + P5'^^ 



a four velocity 



Vl — 



(5) 



(6) 



and a 3- velocity with respect to a zero angular momentum 
observer of 

V = {n - Lo)e'''-'' , (7) 

where and 0" are the two killing vectors associated with 
the time and translational symmetries of the space-time, 
gab is the metric tensor and f2 is the angular velocity. 



Having specified an equation of state for very dense mat- 
ter, the structure of the star is computed by solving four 
components of Einstein's gravitational field equations 



Rab — 8Tr{Tab — -^gabT), 



(8) 



(where Rab is the Ricci tensor and T ^ T^) and the equa- 
tion of hydrostationary equilibrium. 

2.3 Equations of State 

The simplest equation of state one can use to model rel- 
ativistic stars is the relativistic polytropic EOS [p| 



P 



e = pc 



r- 1' 



(9) 
(10) 



where p is the rest mass density, if is a constant and T is 
the polytropic exponent. Instead of F, one often uses the 
polytropic index N , defined through 



1 



1 

N' 



(11) 



For this equation of state, the quantity 
^(r-2)/(r-i)^^i/(r-i)/G has units of length. In 
gravitational units (c = G — 1), one can thus use 
as a fundamental length scale to define dimensionless 
quantities. Equilibrium models are then characterized 
by the polytropic index N and their properties can be 
scaled to different values, using an appropriate value 
for K. For iV < 1.0 (iV > 1.0) one obtains stiff (soft) 
models, while for N = 0.5 — 1.0, one obtains models with 
bulk properties that are comparable to those of observed 
neutron stars. 

Note that the for the above polytropic EOS, the poly- 
tropic index F coincides with the adiabatic index of a 
relativistic isentropic fluid 



F 



e + PdP 

p H' 



(12) 



This is not the case for the polytropic equation of state 
P = Ke^ , that has been used by other authors, which 
satisfies ( p^ ) only in the Newtonian limit. 

The true equation of state that describes the interior 
of compact stars is largely unknown. This results from 
the inability to verify experimentally the different theories 
that describe the strong interactions between baryons and 
the many-body theories of dense matter at densities larger 
than about twice the nuclear density (i.e. at densities 
larger than about 5 x 10^'^gr/cm'^). 

Many different realistic EOSs have been proposed to 
date which all produce neutron stars that satisfy the cur- 
rently available observational constraints (currently, the 
two main constraints are that the EOS must admit non- 
rotating neutron stars with gravitational mass of at least 
IAAMq and allow rotational periods at least as small as 



1.56 ms, see ||, |^). The proposed EOSs are qualitatively 
and quantitatively very different from each other. Some 
are based on relativistic many-body theories while oth- 
ers use nonrelativistic theories with baryon-baryon inter- 
action potentials. A classic collection of early proposed 
EOSs was compiled by Arnett and Bowers [||, while re- 
cent EOSs are described in Salgado et al. 

High density equations of state with pion condensa- 
tion have been proposed by Migdal and Sawyer and 
Scalapino |l^. The possibility of Kaon condensation is 
discussed by Brown and Bethe ||l^ and questioned by 
Pandharipande et al. [Q. Many authors have examined 
the possibility of stars composed of strange quark matter 
and a recent review can be found in [Q . 

The realistic EOSs are supplied in the form of an energy 
density vs. pressure table and intermediate values are in- 
terpolated. This results in some loss of accuracy because 
the usual interpolation methods do not preserve thermo- 



dynamical consistency. Recently however, Swesty 15 de- 



vised a cubic Hermite interpolation scheme that does pre- 
serve thermodynamical consistency and the scheme has 
been shown to indeed produce higher accuracy neutron 
star models in Nozawa et al. fl^ . 

• Going further. A discussion of hybrid stars, that 
have a mixed-phase region of quark and hadronic 
matter, can be found in |l^. A study of the relax- 
ation effect in dissipative relativistic fluid theories is 
presented in Q. 

2.4 Numerical Schemes 

Out of the ten components of the field equations that de- 
scribe the geometry of a rotating relativistic star, only 
four are independent and one has the freedom to choose 
which four components to use. After choosing four field 
equations, there are different methods one can use to 
solve them. First models were obtained by Wilson |19 
and Bonazzola and Schneider pO[ |. Here we will review 
the following methods: Hartlc's slow rotation formalism, 
the Newton-Raphson linearization scheme due to Butter- 
worth and Ipser [pi| , a scheme using Green's functions by 
Komatsu et al. |22|, |2^, a minimal surface scheme due to 
Wu et al. |M, and two spectral methods by Bonazzola 
et al. |2^, |26| Below we give a description about each 
method and its various implementations (codes). 

2.4.1 Hartle 

To 0{il^) the structure of a star changes only by 
quadrupole terms and the equilibrium equations become 



a set of ordinary differential equations. Hartle's |28 
method computes rotating stars in this slow-rotation ap- 
proximation and a review of slowly rotating models has 
been compiled by Datta H. Weber et al. [|o), H] also 
implement Hartle's formalism to explore the rotational 
properties of four new EOSs. 



Weber and Glendenning attempt to improve on 
Hartle's formalism in order to obtain a more accurate 
estimate of the angular velocity at the mass-shedding 
limit but their models show large discrepancies compared 
to corresponding models computed with fully rotating 
schemes Thus, Hartle's formalism cannot be used 
to compute models of rapidly rotating relativistic stars 
with sufficient accuracy. 

2.4.2 Butterworth and Ipser (BI) 

The BTscheme ||2^ solves the four field equations fol- 
lowing a Newton-Raphson like linearization and iteration 
procedure. One starts with a nonrotating model and in- 
creases the angular velocity in small steps, treating a new 
rotating model as a linear perturbation of the previously 
computed rotating model. Each linearized field equation 
is discretized and the resulting linear system is solved. 
The four field equations and the hydrostationary equilib- 
rium equation are solved separately and iteratively until 
convergence is achieved. 

The space is truncated at a finite distance from the 
star and the boundary conditions there are imposed by 
expanding the metric potentials in powers of l/r. Angular 
derivatives are approximated by high-accuracy formulae 
and models with density discontinuities are treated spe- 
cially at the surface. An equilibrium model is specified 
by fixing its rest mass and angular velocity. 

The original BI code was used to construct uniform 
density models and polytropic models ^ . Friedman 
et al. [|4| ^ , extend the BI code to obtain a large number 
of rapidly rotating models based on a variety of realistic 
EOSs. Lattimer et al. used a code which was also 
based on the BI scheme to construct rotating stars using 
recent "exotic" and schematic EOSs, including pion or 
Kaon condensation and self-bound strange quark matter. 

2.4.3 Komatsu, Eriguchi and Hachisu (KEH) 

In the KEH scheme |2^, , the same set of field equations 
as in BI is used, but the three elliptic-type field equations 
are converted into integral equations using appropriate 
Green's functions. The boundary conditions at large dis- 
tance from the star are thus incorporated into the integral 
equations, but the region of integration is truncated at a 
finite distance from the star. The fourth field equation 
is an ordinary first-order differential equation. The field 
equations and the equation of hydrostationary equilib- 
rium are solved iteratively, fixing the maximum energy 
density and the ratio of the polar radius to the equato- 
rial radius, until convergence is achieved. In ^ and 
p7| the original KEH code is used to construct uniformly 
and differentially rotating stars for both polytropic and 
realistic EOSs. 

Cook, Shapiro and Teukolsky (CST) improve on the 
KEH scheme by introducing a new radial variable which 
maps the semi-infinite region [0, oo) to the closed region 



[0,1]. In this way, the region of integration is not trun- 
cated and the model converges to a higher accuracy. De- 
tails of the code are presented in |3q| and polytropic and 
realistic models are computed in [ 39j and Q . 

Stergioulas and Friedman (SF) implement their own 
KEH code following the CST scheme. They improve on 
the accuracy of the code by a special treatment of the 
second order radial derivative that appears in the source 
term of the first-order differential equation for one of the 
metric functions. This derivative was introducing a nu- 
merical error of 1%— 2% in the bulk properties of the most 
rapidly rotating stars computed in the original implemen- 
tation of the KEH scheme. The SF code is presented in 
pl| and in It is available as a public domain code, 
named rns, and can be downloaded from p3| . 

2.4.4 Wu et al. (WMSHR) 

The numerical scheme by Wu et al. implements 
the minimal surface formalism for rotating axisymmetric 
space-times Q, in which Einstein's field equa- 

tions are equivalent to the minimal surface equations in an 
abstract Riemannian potential space with a well-defined 
metric, whose coordinates are the four metric functions 
of the usual stationary, axisymmetric metric. A finite el- 
ement technique is used and the system of equations is 
solved by a Newton-Raphson method. Models based on 
realistic EOSs are presented in The WMSHR code 
has been used to visualize rapidly rotating stars by em- 
bedding diagrams and 4D-ray-tracing pictures (see |48 
for a review). 

2.4.5 Bonazzola et al. (BGSM) 

In the BGSM scheme the field equations are derived 
in the 3-1-1 formulation. All four equations describing 
the gravitational field are of elliptic type. This avoids 
the problem with the second-order radial derivative in 
the source term of the ODE used in BI and KEH. The 
equations are solved using a spectral method, i.e. all func- 
tions are expanded in terms of trigonometric functions in 
both the angular and radial directions and a Fast Fourier 
Transform (EFT) is used to obtain coefficients. Outside 
the star a redefined radial variable is used, which maps 
infinity to a finite distance. 

In the code is used to construct a large number 
of models based on recent EOSs. The accuracy of the 
computed models is estimated using two general relativis- 
tic Virial identities, valid for general asymptotically flat 
space-times, that were discovered by Gourgoulhon and 
Bonazzola |5l| . 

While the field equations used in the BI and KEH 
schemes assume a perfect fiuid, isotropic stress-energy 
tensor, the BGSM formulation makes no assumption 
about the isotropy of Tab- Thus, the BGSM code can 
compute stars with magnetic field, solid crust or solid in- 
terior and it can also be used to construct rotating boson 



stars. 

Since it is based on the 3 + 1 formalism, the BGSM 
code is also suitable for providing high-accuracy, unstable 
equilibrium models as initial data for an axisymmetric 
collapse computation. 

2.4.6 Bonazzola et al. (BGM-98) 

The BGSM spectral method has been improved by Bonaz- 
zola et al. allowing for several domains of integration. 
One of the domain boundaries is chosen to coincide with 
the surface of the star and a regularization procedure is 
introduced for the infinite derivatives at the surface (that 
appear in the density field when stiff equations of state 
are used). This allows models to be computed that are 
free of Gibbs phenomena at the surface. The method is 
also suitable for constructing quasi-stationary models of 
binary neutron stars. 

2.4.7 Direct Comparison of Numerical Codes 

The accuracy of the above numerical codes can be es- 
timated, if one constructs exactly the same models with 
different codes and compares them directly. The first such 
comparison of rapidly rotating models constructed with 
the FIP and SF codes is presented by Stergioulas and 
Friedman in [Q. Rapidly rotating models constructed 
with several EOS's agree to 0.1% — 1.2% in the masses 
and radii and to better than 2% in any other quantity that 
was compared (angular velocity and momentum, central 
values of metric functions etc.). This is a very satisfac- 
tory agreement, considering that the BI code was using 
relatively few grid points, due to limitations of computing 
power at the time of its implementation. 

In it is also shown that a large discrepancy be- 
tween certain rapidly rotating models, constructed with 
the FIP and KEH codes, that was reported by Eriguchi et 
al. was only due to the fact that a different version 
of a tabulated EOS was used in |^ than by FIP. 

Recently, Nozawa et al. |l^ have completed an exten- 
sive direct comparison of the BGSM, SF and the orig- 
inal KEH codes, using a large number of models and 
equations of state. More than twenty different quantities 
for each model are compared and the relative differences 
range from 10^'^ to 10"** or better, for smooth equations 
of state. The agreement is excellent for soft polytropes, 
which shows that all three codes are correct and compute 
the desired models to an accuracy that depends on the 
number of grid-points used to represent the spacetime. 

If one makes the extreme assumption of uniform den- 
sity, the agreement is at the level of lO^'^. In the BGSM 
code this is due to the fact that the spectral expansion in 
terms of trigonometric functions cannot accurately rep- 
resent functions with discontinuous first-order derivatives 
at the surface of the star. In the KEH and SF codes, the 
three-point finite-difference formulae cannot accurately 



represent derivatives across the discontinuous surface of 
the star. 

The accuracy of the three codes is also estimated by 
the use of the two Virial identities due to Gourgoulhon 
and Bonazzola (50[ |l|. Overall, the BGSM and SF codes 
show a better and more consistent agreement than the 
KEH code with BGSM or SF. This is largely due to 
the fact that the KEH code does not integrate over the 
whole spacetime but within a finite region around the 
star, which introduces some error in the computed mod- 
els. 

• Going further. A review of spectral methods in 
general relativity can be found in |52| . A formula- 
tion for nonaxisymmetric, uniformly rotating equi- 
librium configurations in the second post-Newtonian 
approximation is presented in [5^ . 

2.5 Properties of Equilibrium Models 

2.5.1 Bulk Properties of Equilibrium Models 

Neutron star models constructed with various realistic 
EOSs have considerably different bulk properties, due to 
the large uncertainties in the equation of state at high 
densities. Very compressible (soft) EOSs produce models 
with small maximum mass, small radius, and large ro- 
tation rate. On the other hand, less compressible (stiff) 
EOSs produce models with a large maximum mass, large 
radius, and low rotation rate. 

The gravitational mass, equatorial radius and rota- 
tional period of the maximum mass model constructed 
with one of the softest EOSs (EOS B) (l.eSM©, 9.3km, 
0.4ms) are a factor of two smaller than the mass, radius 
and period of the corresponding model constructed by one 
of the stiffest EOSs (EOS L) (3.27Mq, 18.3km, 0.8ms). 
The two models differ by a factor of 5 in central energy 
density and a factor of 8 in the moment of inertia! 

Not all properties of the maximum mass models be- 
tween proposed EOSs differ considerably. For exam- 
ple, most realistic EOSs predict a maximum mass model 
with a ratio of rotational to gravitational energy T/W of 
0.11 ±0.02, a dimensionless angular momentum cJ/GNP 
of 0.64±0.06 and an eccentricity of 0.66±0.04, (l). Hence, 
between the set of realistic EOSs, some properties are 
directly related to the stiffness of the EOS while other 
properties are rather insensitive to stiffness. 

Compared to nonrotating stars, the effect of rotation is 
to increase the equatorial radius of the star and also to 
increase the mass that can be sustained at a given central 
energy density. As a result, the mass of the maximum 
mass rotating model is roughly 15%— 20% higher than the 
mass of the maximum mass nonrotating model, for typical 
realistic EOSs. The corresponding increase in radius is 
30% - 40%. 

The deformed shape of a rapidly rotating star creates a 
distortion, away from spherical symmetry, in its gravita- 
tional field. Far from the star, the distortion is measured 



by the quadrupole- moment tensor Qab- For uniformly ro- 
tating, axisymmetric and equatorially symmetric config- 
urations, one can define a scalar quadrupole moment Q, 
which can be extracted from the asymptotic expansion, 
at large r, of the metric function v. 

Laarakkers and Poisson js^, numerically compute the 
scalar quadrupole moment Q for several equations of 
state, using the rotating neutron star code rns |^ . They 
find that for fixed gravitational mass M, the quadrupole 
moment is given as a simple quadratic fit 



(13) 



where J is the angular momentum of the star and a is 
a dimensionless quantity that depends on the equation 
of state. The above quadratic fit reproduces Q with a 
remarkable accuracy. The quantity a varies between a ^ 
2 for very soft EOSs and a ^ 8 for very stiff EOSs, for 
M = IAMq neutron stars. 

For a given zero-temperature EOS, the uniformly ro- 
tating equilibrium models form a 2-dimensional surface in 
the 3-dimensional space of central energy density, gravi- 
tational mass and angular momentum |4l[| . The surface 
is limited by the nonrotating models (J = 0) and by the 
models rotating at the mass-shedding (Kepler) limit, i.e. 
at the maximum allowed angular velocity so that the star 
does not shed mass at the equator. Cook et al. |40[| 
have shown that the model with maximum angular veloc- 
ity does not coincide with the maximum mass model, but 
is generally very close to it in central density and mass. 
Stergioulas and Friedman ]4l} | show that the maximum 
angular velocity and maximum baryon mass equilibrium 
models are also distinct. The distinction becomes signifi- 
cant in the case where the EOS has a large phase transi- 
tion near the central density of the maximum mass model, 
otherwise the models of maximum mass, baryon mass, an- 
gular velocity and angular momentum can be considered 
to coincide for most purposes. 

2.5.2 An Empirical Formula for the Kepler Ve- 
locity 

In the Newtonian limit the maximum angular veloc- 
ity of uniformly rotating polytropic stars is, ilmax — 
(2/3)3/2(GAf/i?3)i/2 (see 



150 



For relativistic stars. 



the empirical formula |5^, |35| , 56 



0.67. 



IGM„ 



*"max 



(14) 



gives the maximum angular velocity in terms of the mass 
and radius of the maximum mass nonrotating model with 
an accuracy of 5% — 7%, without actually having to con- 
struct rotating models. 

The empirical formula results from universal propor- 
tionality relations that exist between the mass and radius 
of the maximum mass rotating model and those of the 



maximum mass nonrotating model for the same EOS. La- 
sota et al. |5^ find that, for most EOSs, the coefficient in 
the empirical formula is an almost linear function of the 
parameter 

= -5 (15) 

When this relation is taken into account in the empirical 
formula, it reproduces the exact values with a relative 
error of only 1.5%. 



Weber and Glendenning 1 30 , 3^ , try to reproduce ana- 
lytically the empirical formula in the slow rotation ap- 
proximation but the formula they obtain involves the 
mass and radius of the maximum mass rotating config- 



uration, which is different from what is involved in (14) 



2.5.3 The Upper Limit on Mass and Rotation 

The maximum mass and minimum period of rotating rel- 
ativistic stars computed with realistic EOSs from the Ar- 
nett and Bowers collection ||] are about 3.3M0 (EOS L) 
and 0.4ms (EOS B), while 1.4M0 neutron stars, rotat- 
ing at the Kepler limit, have a rotational periods between 
0.53ms (EOS B) and 1.7ms (EOS M) The maxi- 

mum, accurately measured, neutron star mass is currently 
1.44M0, but there are also indications for 2.OM0 neutron 
stars [^8| . The minimum observed pulsar period is 1.56ms 
, which is close to the experimental sensitivity of recent 
pulsar searches (an ongoing experiment is designed to de- 
tect sub- millisecond pulsars, if they exist ||59|). 

In principle, neutron stars with maximum mass or min- 
imum period could exist, if they are born as such in a core 
collapse, or if they accrete the right amount of matter and 
angular momentum during an accretion-induced spin-up 
phase. Such a phase could also follow the creation of an 
I.4M0 neutron star during the accretion induced collapse 
of a white dwarf. 

In reality, only a very small fraction, if any, of neutron 
stars will be close to the maximum mass or minimum pe- 
riod limit. In addition, rapidly rotating nascent neutron 
stars are subject to a nonaxisymmetric instability, which 
lowers their initial rotation rate and neutron stars with 
a strong magnetic field have their rotation rate limited 
by the Kepler velocity at their Alfven radius, where the 
accretion pressure balances the magnetospheric pressure 
§. 

• Going further. A recent review by J. L. Friedman 
on the upper limit on rotation of relativistic stars can 
be found in Jsof . 

2.5.4 The Upper Limit on Mass and Rotation 
Set by Causality 

Current proposed EOSs are reliable only to about twice 
nuclear density and result in very different values for the 
maximum mass and minimum period of neutron stars. 
If one is interested in obtaining upper limits on the mass 



and rotation rate, independent of the proposed EOSs, one 
has to rely on fundamental physical principles. 

Instead of using realistic EOSs, one constructs a set of 
artificial EOSs that satisfy only a minimal set of physi- 
cal constraints, which represent what we know about the 
equation of state of matter with high confidence. One 
then searches among all these EOSs to obtain the one 
that gives the maximum mass or minimum period. The 
minimal set of constraints that have been used in such 
searches are that 

1. the high density EOS matches to the known low den- 
sity EOS at some matching energy density Cm, 

2. the matter at high densities satisfies the causality 
constraint (the speed of sound is less than the speed 
of light). 

In relativistic perfect fluids, the speed of sound is the 
characteristic velocity of the fluid evolution equations and 
the causality constraint translates into the requirement 



models and using the empirical formula (|T^) to estimate 
the minimum period. Koranda, Stergioulas and Friedman 
[ p6[ improve on Glendenning's results by constructing ac- 
curate rapidly rotating models and show that Glenden- 
nings results are accurate to within the accuracy of the 
empirical formula. 

Furthermore, they show that the EOS satisfying the 
minimal set of constraints and yielding the minimum pe- 
riod star consists of a high density region at the causal 
limit, which is matched to the known low density EOS 
through an intermediate constant pressure region (that 
would correspond to a first-order phase transition). Thus, 
the EOS yielding absolute minimum period models is as 
stiff as possible at the central density of the star (to sus- 
tain a large enough mass) and as soft as possible in the 
crust, in order to have the smallest possible radius (and 
rotational period). 

The absolute minimum period of uniformly rotating 
stars is an (almost linear) function of the maximum ob- 
served mass of nonrotating neutron stars 



dp/de < 1. 



(16) 



P™„ = 0.28ms + 0.2(M, 



nonrot. 
max 
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(see e.g. Geroch and Lindblom |61 ). It is assumed that 
the fluid will still behave as a perfect fluid when it is 
perturbed from equilibrium. 

For nonrotating stars, Rhoades and Ruffini showed that 
the EOS that satisfies the above two constraints and 
yields the maximum mass consists of a high density region 
as stiff as possible (i.e. at the causal limit, dp/de = 1), 
that matches directly to the known low density EOS. For 
a chosen matching density e^, they computed a maxi- 
mum mass of 3.2M0. However, this is not the theoreti- 
cally maximum mass of nonrotating neutron stars, as is 
often quoted in the literature. Hartle and Sabbadini 
point out that M^ax is sensitive to the matching energy 



and is rather insensitive to the matching density e„i (the 
above result was computed for a matching number density 
of O.lfm"^). 

In |l66| , it is also shown that an absolute limit on the 
minimum period exists even without requiring that the 
EOS matches to a known low density EOS (this is not true 
for the limit on the maximum mass). Thus, using causal- 
ity as the only constraint on the EOS, Pmin is lowered 
by only 3%, which shows that the currently known part 
of the nuclear EOS plays a negligible role in determin- 
ing the absolute upper limit on the rotation of uniformly 
rotating, gravitationally bound stars. 



density and Hartle |]63| computes Mma 

U4„„ / 3 



as a function of 2.5.5 Spin- Up Prior to Collapse 



. o / 2 X IQi'^gr/cm^ \ . 
M^,„,, = 4.8 ( ^ Mq . 



(17) 



In the case of rotating stars, Friedman and Ipser |6J] 
assume that the absolute maximum mass is obtained by 
the same EOS as in the nonrotating case and compute 
Mmax as a function of matching density, assuming the 
BPS EOS holds at low densities. Stergioulas and Fried- 
man 1^ recompute M™^^ for rotating stars using the 
more recent EPS EOS at low densities, obtaining very 
nearly the same result 



(18) 



where, 2 x 10^''gr/cm'^ is roughly nuclear saturation den- 
sity for the EPS EOS. 

A first estimate of the absolute minimum period of 
uniformly rotating, gravitationally bound stars was com- 
puted by Glendenning p5[ by constructing nonrotating 



Since rotation increases the mass that a neutron star of 
given central density can support, there exist sequences 
of neutron stars with constant baryon number that have 
no nonrotating member. Such sequences are called supra- 
massive as opposed to normal sequences that do have a 
nonrotating member. A nonrotating star can become 
supra-massive by accreting matter and spinning-up to 
large rotation rates; in another scenario, neutron stars 
could be born supramassive after a core collapse. A supra- 
massive star evolves along a sequence of constant baryon 
mass, slowly loosing angular momentum. Eventually, the 
star reaches a point where it becomes unstable to axisym- 
metric perturbations and collapses to a black hole. The 
instability grows on a secular timescale, in the sense that 
it is limited by the time required for viscosity to redis- 
tribute the star's angular momentum. This timescale is 
comparable with the spin- up time following a glitch [Q. 

Cook et al. |3^, Q have discovered that a supra- 
massive star approaching the axisymmetric instability, 
will actually spin-up before collapse, even though it looses 



angular momentum. This, potentially observable, effect 
is independent of the equations of state and it is more pro- 
nounced for rapidly rotating massive stars. In a similar 
phenomenon, normal stars can spin-up by loss of angular 
momentum near the Kepler limit, if the equation of state 
is extremely stiff or extremely soft. 

2.5.6 Rotating Magnetized Neutron Stars 

The presence of a magnetic field was ignored in the models 
of rapidly rotating relativistic stars that were considered 
in the previous sections. The reason is that the observed 
surface dipole magnetic field strength of pulsars ranges 
between B = 10^ G and B = 2x lO^^ q These values of 
B imply a magnetic field energy density that is too small 
compared to the energy density of the fluid, to signifi- 
cantly affect the structure of a neutron star. However, one 
cannot exclude the existence of neutron stars with higher 
magnetic field strengths or the possibility that neutron 
stars are born with much stronger magnetic fields, which 
then decay to the observed values (Of course there are also 
many arguments against magnetic field decay in neutron 
stars 1^). In addition, even though moderate magnetic 
field strengths do not alter the bulk properties of neutron 
stars, they may have an effect on the damping or growth 
rate of various perturbations of an equilibrium star, af- 
fecting its stability. For these reasons, a fully relativistic 
description of magnetized neutron stars is desirable and, 
in fact, Bocquet et al. achieved the first numerical 
computation of such configurations. Here we give a brief 
summary of their work: 

A magnetized relativistic star in equilibrium can be de- 
scribed by the coupled Einstein-Maxwell field equations 
for stationary, axisymmetric rotating objects with inter- 
nal electric currents. The stress-energy tensor includes 
the electromagnetic energy density and is non-isotropic 
(in contrast to the isotropic perfect fluid stress- energy 
tensor) . The equilibrium of the matter is given not only 
by the balance between the gravitational force and the 
pressure gradient, but the Lorentz force due to the elec- 
tric currents also enters the balance. For simplicity, Boc- 
quet et al. consider only poloidal magnetic fields, which 
preserve the circularity of the space-time. Also, they only 
consider stationary configurations, which excludes mag- 
netic dipole moments non-aligned with the rotation axis, 
since in that case the star emits electromagnetic and grav- 
itational waves. The assumption of stationarity implies 
that the fluid is necessarily rigidly rotating (if the mat- 
ter has infinite conductivity) Under these assump- 
tions, the electromagnetic field tensor i^"'' is derived from 
a potential 1-form Aa with only two non-vanishing com- 
ponents, At and A^, which are given by a scalar Poisson 
and a vector Poisson equation respectively. Thus, the 
two equations describing the electromagnetic field are of 
similar type as the four field equations that describe the 
gravitational field. 

The construction of magnetized models with B < 10^^ 



G confirms that magnetic fields of this strength have a 
negligible effect on the structure of the star. However, 
if one increases the strength of the magnetic field above 
10^* G, one observes significant effects, such as a flatten- 
ing of the star. The magnetic field cannot be increased 
indefinitely, but there exists a maximum value of the mag- 
netic field strength, of the order of 10^^ G, for which the 
magnetic field pressure at the center of the star equals 
the fluid pressure. Above this value, the fluid pressure 
decreases more rapidly away from the center along the 
symmetry axis, than the magnetic pressure. Instead of 
pressure, there is tension along the symmetry axis and no 
stationary conflguration can exist. 

The shape of a strongly magnetized star is flattened 
because the Lorentz forces exerted by the E/M field on 
the fluid act as centrifugal forces. A star with a magnetic 
field near the maximum value for stationary configura- 
tions, displays a pinch along the symmetry axis, because 
there, the magnetic pressure exceeds the fluid pressure. 
The maximum fluid density inside the star is not attained 
at the center, but away from it. The presence of a strong 
magnetic field also allows a maximum mass configuration 
with larger Mmax than for the same EOS with no mag- 
netic field and this is in analogy with the increase of Mmax 
induced by rotation. For nonrotating stars, the increase 
in Mraax^ duc to & strong magnetic field, is 13% — 29%, 
depending on the EOS. Following the increase in mass, 
the maximum allowed angular velocity for a given EOS 
also increases in the presence of a magnetic field. 

Bocquet et al. are planning to use their code in the 
study of two types of possible instabilities in magnetized 
neutron stars, i) a pure E/M instability towards another 
electric current/magnetic field distribution of lower en- 
ergy and ii) a nonaxisymmetric instability for rapidly ro- 
tating models, which would be the analog of a Jacobi- 
type transition in non-magnetized stars. In perfect fiuid 
models with a magnetic field, one would also expect a 
CFS-instability driven by electromagnetic waves. 

2.5.7 Rapidly Rotating Proto-Neutron Stars 

Following the gravitational collapse of a massive stellar 
core, a proto-neutron star (PNS) is born. Initially it 
has a large radius of about 100km and a temperature 
of 50-lOOMeV. The PNS may be born with a large rota- 
tional kinetic energy and initially it will be differentially 
rotating. Due to the violent nature of the gravitational 
collapse, the PNS pulsates heavily, emitting significant 
amounts of gravitational radiation. After a few hundred 
pulsational periods, bulk viscosity will damp the pulsa- 
tions significantly. Rapid cooling due to deleptonization 
transforms the PNS to a hot neutron star of T ^ lOMeV 
shortly after its formation. In addition, viscosity reduces 
the differential rotation to a nearly uniform rotation on 
a timescale of seconds and the neutron star becomes 
quasi-stationary. Since the details of the PNS evolution 
determine the exact properties of the resulting cold NSs, 



proto-neutron stars must be modeled realistically in order 
to understand the structure of cold neutron stars. 

Hashimoto et al. |6|] and Goussard et al. Q re- 
cently constructed fully relativistic models of rapidly ro- 
tating, hot proto-neutron stars. The authors use finite- 
temperature EOSs 1?^, to model the interior of PNSs. 
Important parameters, which determine the local state of 
matter but are largely unknown, are the lepton fraction 
Yi and the temperature profile. Hashimoto et al. con- 
sider only the limiting case of zero lepton fraction Yi = 
and classical isothermality, while Goussard et al. consider 
several non-zero values for Yi and two different limiting 
temperature profiles - a constant entropy profile and a 
relativistic isothermal profile. In both |Q and |Q, dif- 
ferential rotation is neglected to a first approximation. 

The construction of numerical models with the above 
assumptions shows that, due to the high temperature and 
the presence of trapped neutrinos, PNSs have a signifi- 
cantly larger radius than cold NSs. These two effects give 
the PNS an extended envelope which, however, contains 
only roughly 0.1% of the total mass of the star. This 
outer layer cools more rapidly than the interior and be- 
comes transparent to neutrinos, while the core of the star 
remains hot and neutrino opaque for a longer time. The 
two regions are separated by the "neutrino sphere" . 

Compared to the T — case, an isothermal EOS with 
temperature of 25MeV has a maximum mass model of 
only slightly larger mass. In contrast, an isentropic EOS 
with a nonzero trapped lepton number features a maxi- 
mum mass model that has a considerably lower mass than 
the corresponding model in the T = case and a stable 
PNS transforms to a stable neutron star. If, however, one 
considers the hypothetical case of a large amplitude phase 
transition which softens the cold EOS (such as a Kaon 
condensate), then Mmax of cold neutron stars is lower 
than M„iax of PNSs and a stable PNS with maximum 
mass will collapse to a black hole after the initial cooling 
period. This scenario of delayed collapse of nascent neu- 
tron stars has been proposed by Brown and Bcthe fl^ 
and investigated by Baumgarte et al. [|72| . 

An analysis of radial stability of PNSs shows that, 
for hot PNSs, the maximum angular velocity star almost 
coincides with the maximum mass star, as is also the case 
for cold EOSs. 

Because of their increased radius, PNSs have a differ- 
ent mass-shedding limit than cold NSs. For an isothermal 
profile, the mass-shedding limit proves to be sensitive to 
the exact location of the neutrino sphere. For the EOSs 
considered in |^ and [|9) PNSs have a maximum angular 
velocity that is considerably less than the maximum an- 
gular velocity allowed by the cold EOS. Stars that have 
nonrotating counterparts (i.e. that belong to a normal 
sequence) contract and speed up while they cool down. 
The final star with maximum rotation is thus closer to the 
mass-shedding limit of cold stars than was the hot PNS 
with maximum rotation. Surprisingly, stars belonging to 
a supra-massive sequence exhibit the opposite behavior. 



If one assumes that a PNS evolves without loosing angu- 
lar momentum or accreting mass, then a cold neutron star 
produced by the cooling of a hot PNS has a smaller an- 
gular velocity than its progenitor. This purely relativistic 
effect was pointed out in |Q and confirmed in ||6^. It 
should be noted here, that a small amount of differential 
rotation significantly affects the mass-shedding limit, al- 
lowing more massive stars to exist than uniform rotation 
allows. Taking differential rotation into account, a more 
recent study by Goussard et al. [D suggests that proto- 
neutron stars created in a gravitational collapse cannot 
spin faster than 1.7 ms. 

3 Oscillations and Stability 

The study of oscillations of relativistic stars has the po- 
tential of yielding important information about both the 
bulk properties and the composition of the interior of the 
star i.e. about the equation of state of matter at very 
high densities, in about the same way that helioseismol- 
ogy is providing us with information about the interior of 
the Sun. In a neutron star - accretion disk system, the 
star-disk interaction can drive oscillations and one of the 
possible explanations for kHz quasi-periodic oscillations 
recently discovered in several X-ray sources are neutron 
star pulsations (for an early proposal that such oscil- 
lations may be observable, see ]75|). 

Neutron star pulsations may be a detectable source of 
gravitational radiation. The pulsations can be excited af- 
ter a core collapse or during the final stages of a neutron 
star binary system coalescence. Rapidly rotating neutron 
stars are unstable to the emission of detectable gravita- 
tional waves for a short time after their formation. The 
identification of gravitational waves produced by a neu- 
tron star can lead to the determination of its mass and 
radius and several such determinations can help recon- 
struct the equation of state of matter at very high energy 
densities | |7^ . 

The oscillations of relativistic stars are actually a 
non-linear phenomenon and their numerical computation 
would require a full 3-D relativistic hydrodynamics code, 
which is not yet available. However, apart from the ini- 
tial oscillations following core collapse, the oscillations of 
an equilibrium star are of small magnitude compared to 
its radius and it will suffice to approximate them as lin- 
ear perturbations. Such perturbations can be described 
in two equivalent ways. In the Lagrangian approach, one 
studies the changes in a given fluid element as it oscillates 
about its equilibrium position. In the Eulerian approach, 
one studies the change in fluid variables at a fixed point 
in space. Both approaches have their strengths and weak- 
nesses. 

In the Newtonian limit, the Lagrangian approach has 
been used to develop variational principles ^ but the 
Eulerian approach proved to be more suitable for numeri- 
cal computations of mode frequencies and eigenfunctions 



|79[ 1^, |l|, |8|, m. Clement H used the Lagrangian 
approach to obtain axisymmetric normal modes of rotat- 
ing stars, while nonaxisymmetric solutions were obtained 
in the Lagrangian approach by Imamura et al. p5[ | and 
in the Eulerian approach by Managan |80| and Ipser and 



Lindblom fSl 



3.1 Quasi-Normal Modes of Oscillation 

The spacetime of a nonrotating star is static and spher- 
ically symmetric. A general linear perturbation can be 
written as a sum of quasi-normal modes that are charac- 
terized by the indices (/, m) of the spherical harmonic Yj™ 
and have angular and time-dependence of the form 



f{r)Yi"'{cose)e" 



(20) 



where Q is a scalar unperturbed quantity, ujp is the angu- 
lar frequency of the mode, as measured by a distant iner- 
tial observer and /(r) represents the radial dependence of 
the perturbation. Normal modes of nonrotating stars are 
degenerate in m and it suffices to study the axisymmetric 
(m — 0) case. 

The perturbation of the metric, Sgab, can be expressed 
in terms of spherical, vector and tensor harmonics. These 
are either of "polar" or "axial" parity. Here, parity is de- 
fined as the change in sign under a combination of reflec- 
tion in the equatorial plane and rotation by tt. A polar 
perturbation has parity (—1)', while an axial perturbation 
has parity (—1)'+^. Because of the spherical background, 
the polar and axial perturbations of a nonrotating star 
are completely decoupled. 

A normal mode solution satisfies the perturbed gravi- 
tational field equations 
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and the perturbation of the conservation of the stress- 
energy tensor 

SiVaT"'') = 0. (22) 

For given {l,m), a solution exists for any value of 
the eigenfrequency ujp and it consists of ingoing- and 
outgoing-wave parts. Outgoing modes are defined by the 
discrete set of eigenfrequencies for which there are no in- 
coming waves at infinity. These are the modes that will 
be excited in various astrophysical situations. 

The main modes of pulsation that are known to exist 
in relativistic stars have been classified as follows (/o and 
To are typical frequencies and damping times of the most 
important modes in the nonrotating limit): 

1. Polar fluid modes 

Are slowly damped modes analogous to the Newto- 
nian fiuid pulsations: 

• /(undamental)-mode: surface mode due to the 
interface between the star and its surroundings 
(/o ~ 2khz, Tn < Isec), 



• p(ressure)-modes: nearly radial (/o > 4kHz, 
To > Is), 

• (/(ravity modes): nearly tangential, only exist 
for finite temperature stars (/o < 500Hz, To > 
5s). 

2. Axial fluid modes 

• r(otation) modes: degenerate at zero-frequency 
for nonrotating stars. In a rotating star, gener- 
ically unstable. Frequencies from zero to kHz, 
growth times inversely proportional to a high 
power of the star's angular velocity. 

3. Polar and axial spacetime modes 

• w(ave)-modes: Analogous to the quasi-normal 
modes of a black hole. High frequency, strongly 
damped modes (/o > 6kHz, To ~ 0.1msec). 



For a more detailed description of various modes see 1 86 



3.2 Effect of Rotation on Quasi-Normal 
Modes 

In a continuous sequence of rotating stars, a quasi-normal 
mode of index I is defined as the mode which, in the 
nonrotating limit, reduces to the quasi- normal mode of 
the same index /. Rotation has several effects on the 
modes of a previously nonrotating star: 

1 . The degeneracy in the index m is removed and a non- 
rotating mode of index I is split into 2^ + 1 different 
(Z, m) modes. 

2. Prograde (m < 0) modes are now different than ret- 
rograde (to > 0) modes. 

3. A rotating "polar" Z-mode consists of a sum of purely 
polar and purely axial terms 



A, 
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that is, rotation couples a polar Z-term to an axial 
I ± 1 term (the coupling to the I + 1 term is, however, 
strongly favored over the coupling to the I — 1 term 
[pi[). Similarly, for a rotating "axial" mode. 
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Frequencies and damping times are shifted. In gen- 
eral, frequencies (in the inertial frame) of prograde 
modes increase, while those of retrograde modes de- 
crease with increasing rate of rotation. 



In rotating stars, quasi-normal modes of oscillation 
have only been studied in the slow-rotation limit, in the 
post-Newtonian and in the Cowling Approximations. The 
solution of the fuUy-relativistic perturbation equations for 
a rapidly rotating star is still a very challenging task and 
only recently they have been solved for zero- frequency 
(neutral) modes i|,|9|. 

• Going further. The equations that describe oscilla- 
tions of the solid crust of a rapidly rotating relativis- 
tic star are derived by Prion in ||9^. The effects of 
superfluid hydrodynamics on the oscillations of neu- 
tron stars are investigated by Lindblom and Mendell 



3.3 Axisymmetric Perturbation 

Along a sequence of nonrotating relativistic stars with 
increasing central energy density, there is always a model 
for which the mass becomes maximum. The maximum 
mass turning point marks the onset of a secular instability 
in the fundamental axisymmetric pulsation mode of the 
star. 

Applying the turning point theorem provided by Sorkin 
]95| , Friedman Ipser and Sorkin |Q show that in the case 
of rotating stars the secular axisymmetric instability sets 
in when the mass becomes maximum along a sequence of 
constant angular momentum. An equivalent criterion is 
provided by Cook et al. Q : the secular axisymmetric 
instability sets in when the angular momentum becomes 
minimum along a sequence of constant rest mass. 

The instability develops on a timescale that is limited 
by the time required for viscosity to redistribute the star's 
angular momentum. This timescale is long compared to 
the dynamical timescale and comparable to the spin-up 
time following a pulsar glitch. When it becomes secularly 
unstable, a star evolves in a quasi-stationary fashion until 
it encounters the dynamical instability and collapses to 
a black hole. Thus, the onset of the secular instability 
to axisymmetric perturbations separates stable neutron 
stars from neutron stars that will collapse to a black hole. 

Goussard et al. extend the stability criterion to 

hot protoneutron stars with nonzero total entropy. In this 
case, the loss of stability is marked by the configuration 
with minimum angular momentum along a sequence of 
both constant rest mass and total entropy. 

In the nonrotating limit, Gondek et al. fz^ compute 
frequencies and eigenfunctions of axisymmetric pulsations 
of hot proto-neutron stars and verify that the secular in- 
stability sets in at the maximum mass turning point, as 
is the case for cold neutron stars. 

• Going further The stabilization of a relativistic 
star, that is marginally stable to axisymmetric per- 
turbations, by an external gravitational field, is dis- 
cussed in i97| . 



3.4 Nonaxisymmetric Perturbations 
3.4.1 Nonrotating Limit 

For a spherical star, it suffices to study the m — ax- 
isymmetric modes of pulsation, since the m ^ modes 
can be obtained by a rotation of the coordinate system. 

Thorne, Campolattaro and Price, in a series of pa- 
pers IOC I, initiated the computation of nonradial 
modes by formulating the problem in the Regge- Wheeler 



(RW) gauge 1 101 1 and numerically computing nonradial 
modes for a number of neutron star models. A variational 
method for obtaining eigenfrequencies and eig enfu nctions 
has been constructed by D etweiler and Ipser |102|. Lind- 
blom and Detweiler [ |l03| explicitly reduced the system 
of equations to four first-order ordinary differential equa- 
tions and obtained more accurate eigenfrequencies and 
damping times for a larger set of neutron star models. 
They later realized that their system of equations is some- 
times singular inside the star and obtained an improved 
set of equations which is free of this singularity |104|. 

Chandrasekhar and Ferrari |]9]]] express the nonradial 
pulsations in terms of a fifth-order system in a diagonal 
gauge, which is independent of fluid variables. They thus 
reformulate the problem in a way analogous to the scat- 
tering of g ravitational waves off a black hole. Ipser and 
Price [ |105| show that in the RW gauge, nonradial pulsa- 
tions can be described by a system of two second-order 
equations, which can also be independent of fluid vari- 
ables. In addition, they find that the diagonal gauge of 
Chandrasekhar and Ferrari has a remaining gauge free- 
dom which, when removed, also leads to a fourth-order 



system of equations [106|. 

In order to locate purely outgoing- wave modes, one has 
to be able to distinguish the outgoing-wave part from the 
ingoing- wave part at infinity. In the Thorne et al. and 
Lindblom and Detweiler schemes, this is achieved using 
analytic approximations of the solution at infinity. 

VF-modes pose a more challenging numerical problem 
because they are strongly damped and the techniques 
used for / and p modes fail to distinguish the outgoing- 
wave part, but Andersson, Kokkotas and Schutz [107|, 
successfully combine a redefinition of variables with a 
complex-coordinate integration method, obtaining highly 
accurate complex frequencies for w modes. In this 
method, the ingoing and outgoing solutions are separated 
by numerically calculating their analytic continuations to 
a place in the complex-coordinate place, where they have 
comparable amplitudes. Since this approach is purely nu- 
merical, it could prove to be suitable for the computation 
of quasi-normal modes in rotating stars, where analytic 
solutions at infinity are not available. 

The non-availability of asymptotic solutions at infinity 
in the case of rotating stars is one of the major diffi- 
culties for computing outgoing modes in rapidly rotating 
relativistic stars. A new development that may help to 
overcome this problem, at least to an acceptable approx- 



imation, is presented in |10S| by Lindblom, Mendell and 
Ipser. 

The authors obtain approximate near-zone bound- 
ary conditions for the outgoing modes that replace the 
outgoing-wave condition at infinity and that enable one 
to compute the eigenfrequencies with very satisfactory 
accuracy. First, the pulsation equations of polar modes 
in the Regge- Wheeler gauge are reformulated as a set 
of two second-order radial equations for two potentials - 
one corresponding to fluid perturbations and the other to 
the perturbations of the spacetime. The equation for the 
space-time perturbation reduces to a scalar wave equation 
at infinity and to Laplace's equation for zero- frequency so- 
lutions. From these, an approximate boundary condition 
for outgoing modes is constructed and imposed in the 
near zone of the star (in fact on its surface) instead at 
infinity. For polytropic models, the near-zone boundary 
condition yields /-mode eigenfrequencies with real parts 
accurate to 0.01% — 0.1% and imaginary parts with ac- 
curacy at the 10% — 20% level, for the most relativistic 
stars. If the near zone boundary condition can be applied 
to the oscillations of rapidly rotating stars, the resulting 
frequencies and damping times should have comparable 
accuracy. 



3.4.2 Slow Rotation Approximation 

The slow rotation approximation has proven to be useful 
for obtaining a first estimate of the effect of rotation on 
the pulsations of relativistic stars. To lowest order in 
rotation, a polar ^-mode of an initially nonrotating star 
couples to an axial I ± 1 mode in the presence of rotation. 
Conversely, an axial l-mode couples to a polar / ± 1 mode 
|9ll. 

The equations of nonaxisymmetric perturbations in the 
slow-rotation limit and in the Regge- Wheeler gauge are 
derived by Kojima in ]109|, 110|, where the complex fre- 



quencies a = afj+iai for the I = m modes of various poly- 
tropes are computed. For countcrrotating modes, both 
an and tr/ decrease, tending to zero, as the rotation rate 
increases (when a passes through zero, the star becomes 
unstable to the CFS-instability). Extrapolating ajj and 
CTj to higher rotation rates, Kojima finds a large discrep- 
ancy between the points where an and aj go through 
zero. This shows that the slow rotation formalism cannot 
accurately determine the onset of the CFS-instability of 
polar modes in rapidly rotating neutron stars. 



In I 111 I, it is shown that, for slowly rotating stars, the 
coupling between polar and axial modes affects the fre- 
quency of pulsation only to second order in rotation, so 
that, in the slow rotation approximation, to 0(fi), the 
coupling can be neglected when computing frequencies. 

The slow rotation approximation has also been used 
recently in the study of the r-mode instability [139, 112 



3.4.3 Post-Newtonian Approximation 

A first step towards the solution of the perturbation equa- 
tions i n fu l l relativit y has been taken by Cutler and Lind- 
blom |lT|, 
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115 ] who obtain frequencies for the / 



/-modes in rotating stars in the first post-Newtonian 
(1-PN) approximation. The perturbation equations are 
derived in the post-Newtonian formalism of Gunnarsen 
[116 1, i.e. the equations are separated into equations of 
consistent order in 1/c. 

Cutler and Lindblom show that in this scheme, the per- 
turbation of the 1-PN correction of the four- velocity of the 
fluid can be obtained analytically in terms of other vari- 
ables, similarly to what is done for the perturbation in the 
four-velocity in the Newtonian Ipser-Managan scheme. 
The perturbation in the 1-PN corrections are obtained 
by solving an eigenvalue problem, which consists of three 
second order equations, with the 1-PN correction to the 
eigenfrequency of a mode, Ato, as the eigenvalue. 

Cutler and Lindblom obtain a formula that yields Auj 
if one knows the 1-PN stationary solution and the solu- 
tion to the Newtonian perturbation equations. Thus, the 
frequency of a mode in the 1-PN approximation can be 
obtained without actually solving the 1-PN perturbation 
equations numerically. The 1-PN code was checked in the 
nonrotating limit and it was found to reproduce the exact 
general relativistic frequencies for stars with M/R = 0.2 
obeying an = 1 polytropic EOS with an accuracy of 
3% - 8%. 

Along a sequence of rotating stars, the frequency of a 
mode is commonly described by the ratio of the frequency 
of the mode in the comoving frame to the frequency of the 
mode in the nonrotating limit. For an = 1 polytrope 
and for AI/R = 0.2, this frequency ratio is reduced by 
as much as 12% in the 1-PN approximation compared to 
its Newtonian counterpart (for the fundamental I = m 
modes) which is representative of the effect that general 
relativity has on the frequency of quasi-normal modes in 
rotating stars. 

3.4.4 Cowling Approximation 



In several situations, the frequency of pulsations in rel- 
ativistic stars can be estimated even if one completely 
neglects the perturbation in the gravitational field, i.e. 
if one sets Sgab = in the perturbation equations [117|. 
In this approximation, the pulsations are described only 
by the perturbation in the fluid variables and the scheme 



works quite well for /, p and r-modes [118|. A different 
version of the Cowling approximation, in which 6gtr is 
kept nonzero in t he perturbation equations, works better 
for g- modes [119|. 



Yoshida and Kojima [120|, examine the accuracy of 
the relativistic Cowling approximation in slowly rotat- 
ing stars. The first-order correction to the frequency of 
a mode depends only on the eigenfrequency and eigen- 
functions of the mode in the absence of rotation and on 



the angular velocity of the star. The eigenfrequencies of 
/, pi and p2 modes for slowly rotating stars with M/R 
between 0.05 and 0.2 are computed (assuming polytropic 
EOSs with = 1 and N — 1.5 and compared to their 
counterparts in the slow- rotation approximation. 

For the I = 2 /-mode, the relative error in the eigenfre- 
quency because of the Cowling approximation is 30% for 
less relativistic stars {M/R = 0.05) and about 15% for 
stars with M/R = 0.2 and the error decreases for higher 
Z-modes. For the pi and p2 modes the relative error is 
similar in magnitude but it is smaller for less relativistic 
stars. Also, for p-modes, the Cowling approximation be- 
comes more accurate for increasing radial mode number. 

As an application, Yoshida and Eriguchi |121| use the 
Cowling approximation to estimate the onset of the CFS 
instability in rapidly rotating relativistic stars. 

3.5 Nonaxisymmetric Instabilities 
3.5.1 Introduction 

Rotating cold neutron stars, detected as pulsars, have 
a remarkably stable rotation period. But, at birth, or 
during accretion, rapidly rotating neutron stars can be 
subject to various nonaxisymmetric instabilities, which 
will affect the evolution of their rotation rate 

If a protoneutron star has a sufficiently high rotation 
rate (larger than T/W ~ 0.27 for uniformly rotating, con- 
stant density Maclaurin spheroids), it will be subject to a 
dynamical instability driven by hydrodynamics and grav- 
ity. Through the I — 2 mode, the instability will deform 
the star into a bar shape. This highly nonaxisymmetric 
configuration will emit strong gravitational waves with 
frequencies in the kHz regime. The development of the 
instability and the resulting waveform have been com- 
puted numerically in the context of New tonian gravity 
and hydrodynamics by Houser et al. |122|. 

At lower rotation rates, the star can become unstable 
to secular nonaxisymmetric instabilities, driven by grav- 
itational radiation or viscosity. Gravitational radiation 
drives a nonaxisymmetric instability when a mode that is 
retrograde with respect to the star appears as prograde 
to a distant observer, via the Chandrasekhar-Friedman- 
Schutz (CFS) mechanism [123, ItS]: A mode that is ret- 



rograde in the corotating frame has negative angular mo- 
mentum, because the perturbed star has less angular mo- 
mentum than the unperturbed one. If, to a distant ob- 
server, the mode appears prograde, it removes positive 
angular momentum from the star and thus the angular 
momentum of the mode becomes increasingly negative. 

The instability evolves on a secular timescale, during 
which the star loses angular momentum via the emitted 
gravitational waves. When the star rotates slow enough, 
the mode becomes stable and the instability proceeds on 
the longer timescale of the next unstable mode, unless it 
is suppressed by viscosity. 

Neglecting viscosity, the CFS-instability is generic in 



rotating stars for both polar and axial modes. For po- 
lar modes, the instability occurs only above some critical 
angular velocity, where the frequency of the mode goes 
through zero in the inertial frame. The critical angular 
velocity is smaller for increasing mode number I. Thus, 
there will always be a high enough mode number I, for 
which a slowly rotating star will be unstable. Axial modes 
are generically unstable in all rotating stars, since the 
mode has zero frequency in the inertial frame when the 
star is nonrotating |139, 140 1. 

The shear and bulk viscosity of neutron star matter 
is able to suppress the growth of the CFS-instability ex- 
cept when the star passes through a certain temperature 
window. In Newtonian gravity, it appears that the po- 
lar mode CFS-instability can occur only in nascent neu- 
tron stars that rotate close to the mass-shedding limit 
[|2|, |8^, 124, 125, 126 1, but the determination of neutral 
/-modes in full relativity |9^ shows that relativity en- 
hances the instability, allowing it to occur in stars with 
smaller rotation rates than previously thought. 

• Going further. A new numerical method for the 
analysis of the ergoregion instability in relativistic 
stars, which may also be used for the analysis of non- 
axisymmetric instabilities, is presented by Yoshida 



abd Eriguchi in [128 1 



3.5.2 CFS-Instability of Polar Modes 

The existence of the CFS-instability in rota ting stars was 
first demonstrated by Chandrasekhar [123 in the case of 
the I = 2 mode in uniformly rotating, constant density 
Maclaurin spheroids. Friedman and Schutz [|7^, show 
that this instability also appears in compressible stars and 
that all rotating self-gravitating perfect fluid configura- 
tions are generically unstable to the emission of gravita- 
tional waves. In addition, they find that a nonaxisymmet- 
ric mode becomes unstable when its frequency vanishes in 
the inertial frame. Thus, zero-frequency outgoing-modes 
in rotating stars are neutral (marginally stable). 

In the Newtonian limit, neutral modes have been deter- 
mined for several polytropic EOSs g |8^, |l[ |l2|. The 
instability first sets in through I = m modes. Modes with 
larger I become unstable at lower rotation rates but vis- 
cosity limits the interesting ones to / < 5. For an = 1 
polytrope, the critical values of T/W for the I = 3,4 and 
5 modes are 0.079, 0.058 and 0.045 respectively and these 
values become smaller for softer polytropes. 

The I = m = 2 "bar" mode behaves considerably dif- 
ferent than the other modes. It's critical T/W ratio is 
0.14 and it is almost independent of the polytropic in- 
dex. Since soft EOSs cannot produce models with high 
T/W values, the bar mode instability appears only for 
stiff Newtonian polytropes of < O.i 



[129, 130 



In 

addition, the viscosity driven bar mode appears at the 
same critical T/W ratio as the bar mode driven by grav- 
itational radiation (we will see later that this is no longer 
true in general relativity). 



The post-Newtonian computation of neutral modes by 



Cutler and Lindblom |114, 115| has shown that general 
relativity tends to strengthen the CFS-instability. Com- 
pared to their Newtonian counterparts, critical angular 
velocity ratios ilc/^o (where fta = (3Mo/4i?j^)^/^ and 
Mq, Rq are the mass and radius of the nonrotating star 
in the sequence), are lowered by as much as 10% for stars 
obeying the iV = 1 polytropic EOS (for which the in- 
stability occurs only for I = m > 3 modes in the post- 
Newtonian approximation) . 

In full general relativity, neutral modes have been de- 
termined for polytropic EOSs of iV > 1.0 by Stergioulas 
and Friedman p^ , p2| , using a new numerical scheme. 
The scheme completes the Eulerian formalism developed 
by Ipser and Lindblom in the Cowling approximation 
(where Sgab was neglected) |124], by finding an appro- 



priate gauge in which the time-independent perturbation 
equations can be solved numerically for Sgab- Because 
linear perturbations have a gauge freedom, four out of 
ten components of Sgat are fixed by the choice of gauge. 
In the Ipser and Lindblom scheme, the perturbed Eu- 
ler equations are solved analytically. A complete neutral 
mode solution of the perturbation equations is then deter- 
mined by setting the frequency in the inertial frame equal 
to zero and solving six perturbed field equations for Sgab 
and the perturbed equation of energy conservation for a 
scalar function SU. 

The six perturbed field equations in the gauge of Ster- 
gioulas and Friedman are of different types. Three are 
second order ODEs, two are elliptic and the other one is 
parabolic. Their solutions vanish at the center, at infinity 
and on the axis of symmetry, while they are either odd 
or even under reflection in the equatorial plane. The six 
equations, although of different type, are solved simul- 
taneously on a two-dimensional grid, which extends to 
infinity by a redefinition of the radial variable. Solutions 
of the perturbed field equations are obtained for a set of 
trial functions SUi. In the Newtonian limit, it was found 
that the real eigenfunctions can be expanded accurately 
in terms of these trial functions ||8l| . 

The remaining equation to be satisfied, the per- 
turbed energy conservation equation, can be represented 
schematically as a linear operator L on the eigenfunction 
SU. Defining an inner product < SUj\L\SUi >, for the 
set of trial functions, the perturbed energy conservation 
equation is satisfied, when 



det < SUj\L\SU^ >= 0. 



(25) 



Using this criterion, one starts with slowly rotating con- 
figurations and increases the angular velocity of the star 
until (|2^) is satisfied and a complete neutral mode solu- 
tion is obtained. 

The determination of neutral modes for N — 1.0, 1.5 
and 2.0 relativistic polytropes shows that relativity sig- 
nificantly strengthens the instability (which was already 
indicated in the post-Newtonian approximation). For the 



= 1.0 polytrope, the critical angular velocity ratio 
ilc/^K, where flK is the angular velocity at the mass- 
shedding limit at same central energy density, drops by 
as much as 15% for the most relativistic configuration. 
This is a large decrease compared to the Newtonian val- 
ues, which significantly moves the onset of the instability 
away from the mass-shedding limit and which strengthens 
it with respect to the damping effect of viscosity. 

A surprising result, which was not detected in the post- 
Newtonian approximation, is that the I = m ~ 2 bar 
mode is unstable for relativistic polytropes of index A^ = 
1.0. The classical Newtonian result for the onset of the 
bar mode instability {Ncrit < 0.808) is replaced by 



(26) 



Nc„t < 1.3, 



in general relativity. 

Also, in relativistic stars, the onset of the gravitational 
radiation driven bar mode is different from the onset of 
the viscosity driven bar mode. While in the Newtonian 
limit the two bar modes occur at the same critical rotation 
ratio relativity strengthens the gravitational radia- 
tion instability, allowing softer configurations to become 
unstable and suppresses the viscosity driven instability 
allowing it to occur only for very stiff EOSs [131|. 

An independent determination of the onset of the CFS- 
instability in the relativistic Cowling approximation by 
Yoshida and Eriguchi [[l2l| agrees qualitatively with the 



conclusions in |92 

Morsink, Stergioulas and Blattning [132| extend the 
method presented in |Q to a wide range of realistic equa- 
tions of state (which usually have a stiff high density re- 
gion, corresponding to polytropes of index N = 0.5 — 0.7) 
and find that the I ^ m = 2 bar mode becomes unstable 
for stars with gravitational mass as low as 1.0 — 1.2Mo. 
For 1.4Mo neutron stars, the mode becomes unstable at 
80% — 95% of the maximum allowed rotation rate. For a 
wide range of equations of state, the I = m = 2 /-mode 
becomes unstable at a ratio of rotational to gravitational 
energies T/W ~ 0.08 for 1.4Mo stars and T/W - 0.06 
for maximum mass stars. This is to be contrasted with 
the Newtonian value of T/W ^ 0.14. The empirical for- 
mula 

M 

(T/W^)2 = 0.115 - 0.048—^, (27) 

where M^^^ is the maximum mass for a spherical star 
allowed by a given equation of state, gives the critical 
value of T/W for the bar /—mode instability, with an 
accuracy of 4%— 6%, independent of the equation of state. 

Conservation of angular momentum and the inferred 
initial period (assuming magnetic braking) of 6 — 9ms 
for the X-ray pulsar in the supernova remnant N157B 
[133 1, suggests that a fraction of neutron stars may be 
born with very large rotational energies. The /-mode bar 
CFS-instability thus appears as a promising source for the 
planned gravitational wave detectors |134|. It could also 
play a major role in the rotational evolution, through the 



emission of gravitational waves, of merged binary neutron 
stars, if their post-merger angular momentum exceeds the 
maximum allowed to form a Kerr black hole |135|. 



3.5.3 CFS-Instability of Axial Modes 

In nonrotating stars, axial fluid modes are degenerate at 
zero-frequency but in rotating stars they have nonzero 
frequency and are called r-modes in the Newtonian limit 
136, 137|. To 0{fl), their frequency in the inertial frame 



IS 



-mQ 1 



and modes with different radial eigenfunctions can be 
computed at order Q'^ 



112, h3§. According to (M), 



modes with to > are prograde {uji < 0) with respect to 
a distant observer but retrograde {tUr — uJi + mfl > 0) 
in the comoving frame for all values of the angular ve- 
locity. Thus, r-modes in relativistic stars are generically 
unstable to the emission of gravitational waves via the 
CFS-instability, which was first discovered by Andersson 
139 , in the case of slowly rotating, relativistic stars. This 



result is confirmed analytically by Friedman and Morsink 
140 , who show that the canonical energy of the modes 



is negative. 

Two independent computations in the Newtonian 
Cowling approximation, by Andersson, Kokkotas and 
Schutz [[l27| and Lindblom, Owen and Morsink |l4l| show 
that viscosity is not able to damp the r-mode instabil- 
ity in rotating stars. In a temperature window of 10^ K 
< T < 10^° K, the growth time of the ^ = to = 2 mode 
becomes shorter than the shear or bulk viscosity damping 
time at a critical rotation rate that is roughly one tenth 
the maximum allowed angular velocity of uniformly ro- 
tating stars. The gravitational radiation is dominated 
by the current quadrupole term. These results suggest 
that a rapidly rotating proto-neutron star will spin down 
to Crab-like rotation rates within one year of its birth, 
because of the r-mode instability. The current uncertain- 
ties in the viscosity and superfiuid mutual friction damp- 
ing times make this scenario also consistent with some- 
what higher initial spins, like the suggested initial spin 
of 6 — 9 ms o f the X-ray pulsar in the supernova remnant 
N157B [133|. Millisecond pulsars with periods less than 
~ 5ms can then only form after the accretion-induced 
spin-up of old pulsars and not in the accretion-induced 
collapse of a white dwarf. 

The precise limit on the angular velocity of newly-born 
neutron stars will depend on several factors, such as the 
strength of the bulk viscosity, the cooling process, the 
superfiuid mutual friction etc. In the uniform density 
approximation, the r-mode instability can be studied an- 
alytically to 0(r2^) in the angular velocity of the star 
and the resulting expressions for the timescales, given in 
Kokkotas and Stergioulas |142], can be used to study the 
effect of such factors on the instability. In |142| it is also 
shown that the minimum critical angular velocity for the 



onset of the r-mode instability is rather insensitive to the 
choice of equation of state. 

A first study on the issue of detectability of gravita- 
tional waves from the r-mode instability, is presented in 



[144| (see se ctio n 3.5.6), while Andersson, Kokkotas and 
Stergioulas |l43|| study the relevance of the r-mode insta- 
bility in limiting the spin of recycled millisecond pulsars. 

3.5.4 Effect of Viscosity on CFS-Instability 

In the previous sections, we have discussed the growth of 
the CFS-instability driven by gravitational radiation in 
an otherwise nondissipative star. The effect of neutron 
star matter viscosity on the dynamical evolution of non- 
axisymmetric perturbations can be considered separately, 
when the timescale of the viscosity is much longer than 
the oscillation timescale. If tqj^ is the computed growth 
rate of the instability in the absence of viscosity and Tj , 
Tb are the timescales of shear and bulk viscosity, then the 
total timescale of the perturbation is 

i^^ + l-fl. (29) 

Since tgr < and Tb, > 0, a mode will grow only if tgr 
is shorter than the viscous timescales, so that 1/r < 0. 

The shear and bulk viscosity are sensitive to several 
factors and we give here a summary of what is known to 
date from Newtonian and post-Newtonian computations: 

• Shear viscosity 

In normal neutron star matter, shear viscosity is 
dominated by neutron-neutron scattering with a 
temperature dependence of T"^ [145| and compu- 
tations in the Newtonian limit and post-Newtonian 
approximation show that the CFS-instability is sup- 

K II H, |l2| 



pressed for T < 10^ K - 10^ K [|2|, |3|, 

If neutrons become a superfiuid below a transition 
temperature Tg, then mutual friction, which is caused 
by the scattering of electrons off the cores of neutron 
vortices can completely suppress the instability for 
T < Tg. The superffuid transition temperature de- 
pends on the theoretical model for superffuidity and 
lies in the range 10^ K - 6 x 10^ K |l4^ . 



115| 



• Bulk Viscosity 

In a pulsating fiuid that undergoes compression and 
expansion, the weak interaction requires a relatively 
long time to re-establish equilibrium. This creates 
a phase lag between density and pressure perturba- 
tions, which results in a large bulk viscosity [147|. 
The bulk viscosity due to this effect can suppress the 
CFS-instability only for temperatures for which mat- 
ter has become transparent to neutrinos. [ 134 1 , [ 148 1 . 
It has been proposed that for T > 5 x lO^K, matter 
will be opaque to neutrinos and the neutrino phase 
space could be blocked (|134| see also [148|). In this 
case, bulk viscosity will be too weak to suppress the 
instability, but a more detailed study is needed. 



In the neutrino transparent regime, the effect of bulk vis- 
cosity on the instabihty depends crucially on the proton 
fraction Xp. If Xp is lower than a critical value (~ g), only 
modified URCA processes are allowed and bulk viscosity 
limits, but does not suppress completely, the instability 
82, 83, 125 1 . For most modern EOSs, however, the proton 
fract ion is larger than ~ | at sufficiently high densities 



nonaxisymmetric terms are taken into account. 



151 , allowing direct URCA processes to take place. In 
this case, depending on the EOS and the central density 
of the star, the bulk viscosity could almost completely 
suppress the CFS-instability in the neutrino transparent 
regime [152| (but it will probably still not affect it for 
temperatures T > 5 x 10® K). 

In conclusion, the available Newtonian computations 
indicate that the CFS-instability in /—modes is effective 
in nascent neutron stars for temperatures between lO^K 
and lO^^K and possibly also above lO^^K if the star is 
opaque to neutrinos and the bulk viscosity is weak. If di- 
rect URCA reactions do participate in the cooling process, 
it appears that the instability can grow only for temper- 
atures for which the star is opaque to neutrinos. Since 
the neutral mode computations in fully relativistic stars 
show that relativity strengthens the instability, the above 
conclusion should also hold in relativistic stars. 

The uncertainties regarding the effect of viscosity on 
the CFS-instability in realistic neutron stars will be 
greatly reduced by the construction of mode eigenfunc- 
tions for fully relativistic, rotating stars. 

3.5.5 Viscosity-Driven Instability 

A different type of nonaxisymmetric instability in rotat- 
ing stars is that driven by viscosity, which breaks the 
circulation of the fluid [153, 129|. The instability is sup- 
pressed by gravitational radiation, so it can act only 
in cold neutron stars that become rapidly rotating by 
accretion-induced spin-up. The instability sets in when 
the frequency of an Z = — m mode goes through zero in 
the rotating frame. In contrast to the CFS-instability, 
the viscosity-driven instability is not generic in rotating 
stars. The m = 2 mode becomes unstable at a high rota- 
tion rate for very stiff stars and higher m-modes become 
unstable at larger rotation rates. 

In Newtonian polytropes, the instability occurs only for 
stiff polytropes of index N < 0.808 |l29[ |l3§. For rela- 
tivistic models, the situation for the instability becomes 
worse, since relativistic effects tend to suppress the viscos- 
ity instability (while they strengthen the CFS-instability). 
According to recent results by Bonazzola et al. |l3l| , for 
the most relativistic stars, the viscosity driven bar mode 
can become unstable only if < 0.55. For IAMq stars, 
the instability is present for N < 0.67. 

These results are based on an approximate computa- 
tion of the instability in which one perturbs an axisym- 
metric and stationary configuration and studies its evolu- 
tion by constructing a series of triaxial quasi-equilibrium 
configurations. During the evolution only the dominant 



The method presented in |131 is an improvement (tak- 
ing into account nonaxisymmetric terms of higher order) 
of an earlier method by the same authors |148|. Although 
the method is approximate, its results indicate that the 
viscosity-driven instability is likely to be absent in most 
relativistic stars, unless the EOS turns out to be unex- 
pectedly stiff. 

An investigation of the viscosity-driven bar mode in- 
stability, using incompressible, uniformly rotating triax- 
ial ellipsoids in the post-Newtonian approximation, by 
Shapiro and Zane |14£], also finds that the relativistic 
effects weaken the instability. 



3.5.6 Gravitational 
Instability 



Radiation from CFS- 



The CFS-instability can limit the maximum angular ve- 
locity of nascent neutron stars, but it is also a mecha- 
nism for the generation of gravitational waves that could 
be strong enough to be detected by the planned gravita- 
tional wave detectors. 



Lai and Shapiro [134| have studied the development of 
the /-mode instability using Newtonian ellipsoidal rotat- 
ing models |154, |155{ . They consider the case where a 
rapidly rotating neutron star is created in a core collapse. 
After a brief dynamical phase, the protoneutron star be- 
comes axisymmetric but secularly unstable. The instabil- 
ity deforms the star into a nonaxisymmetric configuration 
via the I = 2 bar mode. Since the star looses angular mo- 
mentum via the emission of gravitational waves, it spins- 
down until it becomes secularly stable. 

The frequency of the waves sweeps downward from a 
few hundred Hz to zero, passing through LIGO's ideal 
sensitivity band. A rough estimate of the wave amplitude 
shows that, at ~ lOOHz, the gravitational waves from the 
CFS-instability could be detected out to the distance of 
140Mpc by the advanced LIGO detector. This result is 
very promising, especially since for relativistic stars the 
instability will be stronger than the present Newtonian 
estimate. 

The recently discovered CFS-instability in r-modes, is 
also an important source of gravitational waves. Owen 
et al. 1 144 model the development of the instability and 
the evolution of the neutron star during its spin-down 
phase. The evolution suggests that a neutron star formed 
in the Virgo cluster could be detected by the advanced 
LIGO and VIRGO gravitational wave detectors, with an 
amplitude signal-to-noise ratio that could be as large as 
about 8, if near-optimal data analysis techniques are de- 
veloped. Assuming a substantial fraction of neutron stars 
are born with spin frequencies near their maximum val- 
ues, the stochastic background of gravitational waves pro- 
duced by the r-mode radiation from neutron star forma- 
tion throughout the universe is shown to have an energy 
density of about 10^® of the cosmological closure density, 
in the range 20 Hz to 1 kHz. This radiation is potentially 



detectable by the advanced LIGO as well. 

In newly born stars or in the post-merger objects in 
binary neutron star mergers, rotating close to the Kepler 
limit, both the / and r modes will be unstable. Rcla- 
tivistic computations of growth times in rapidly rotating 
stars or even nonlinear evolutions, are needed to deter- 
mine which mode will be strongest. 

• Going further. The possible ways for neutron stars 
to emit gravitational waves and on their detectability 
are reviewed by Bonaz zola , Gourgou l hon, Flanagan , 
Thorne and Schutz in (l5|, |l57l [Si |l59l pol pi). 
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